Vibrational excitation of diatomic molecular ions in strong-field ionization of diatomic 

molecules 
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A model based on the strong-field and Born-Oppenheimer approximations qualitatively describes 
the distribution over vibrational states formed in a diatomic molecular ion following ionization of 
the neutral molecule by intense laser pulses. Good agreement is found with a recent experiment 
[X. Urbain et al, Phys. Rev. Lett. 92, 163004 (2004)]. In particular, the observed deviation 
from a Franck-Condon-like distribution is reproduced. Additionally, we demonstrate control of the 
vibrational distribution by a variation of the peak intensity or a change of frequency of the laser 
pulse. 

PACS numbers; 33.80.Rv,33.80.Eh,82.50.Hp 



Intense investigations during the last decades have led 
to quite a detailed understanding of the interaction be- 
tween strong laser fields and atoms (see Q for a review). 
This includes, e.g., the process of high-harmonic genera- 
tion which is by now a valuable source for the production 
of coherent ultraviolet light. For molecules, however, the 
extra degrees of freedom introduced by the presence of 
more than one nucleus lead to a much more involved pic- 
ture, and the description of the strong-field driving of 
such systems far from equilibrium is still a challenge in 
theoretical physics (for reviews see 

The main topic of this work is to investigate the par- 
titioning of energy among the electronic and nuclear de- 
grees of freedom when ionizing a diatomic molecule by 
a strong laser field. This is a topic of much current 
interest: until recently it was assumed that the distri- 
bution over vibrational states formed in the molecular 
ion following ionization of the neutral molecule follows a 
Franck-Condon (FC) distribution, and this assumption 
is still frequently used (see jH and references therein). A 
recent experiment, however, reported a non-FC distribu- 
tion and suggested that the use of the FC principle is 
inaccurate because the rate for tunneling ionization in- 
creases sharply with internuclear distance 0. We will 
show that the rapid response of the electron to the laser 
field in fact ensures the applicability of the FC principle 
for the nuclear motion at each instant of time during the 
pulse. It is instead the variation of the ionization and 
excitation rate with intensity in the focal volume that 
leads to a departure from a conventional FC-like distri- 
bution. The theory is very versatile and readily adapted 
to a wide range of diatomic molecules. 

We note that as the strong field response of the molec- 
ular ion is very dependent on the distribution over vibra- 
tional levels, it is desirable to be able to control the latter. 
We show by explicit examples in H2, O2 and N2 that a 
simple change of the peak intensity of a 45 fs Gaussian 
pulse leads to a significant degree of control over the fi- 
nal vibrational distribution, and hence holds the promise 
for the production of a target of interest for state-specific 



studies. 

The energy absorbed by a molecule in the field is dis- 
tributed among the electrons and nuclei. Ionization of 
molecules is accordingly accompanied by a vibrational ex- 
citation of the nuclear motion in the ion formed. To illus- 
trate the accuracy of our approach, we first consider H2 
for which strong-field experiments are available ^ IH 1^ - 
For II2 both photoelectron spectroscopy with UV light 
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sources .7] and ionization by 100 eV electron impact 
resulted in vibrational distributions of in accordance 
with the FC principle, i.e., governed by the projection of 
the initial vibrational ground state wave function v' ^ 
onto the set of vibrational states of the molecular ion, ly. 
In recent experiments 0, 0, 01 > the vibrational distribu- 
tion of H J was measured after ionization of II2 by intense 
laser fields. These experiments revealed very different vi- 
brational distributions than the experiments mentioned 
above. The laser-induced ionization leads to a narrower 
distribution with almost no population in the final vi- 
brational levels above 1/ = A, and the actual shape of 
the distribution depends on the peak laser intensity and 
wavelength. At low intensities, the maximum population 
is observed in the v = state while the FC principle pre- 
dicts a maximum population in the ly = 2 state. These re- 
sults can be partly reproduced by a tunneling model with 
an electronic binding energy which depends on the inter- 
nuclear distance ^Eil^' However, the experimental con- 
ditions in some of the experiments ^ do not correspond 
to a pure tunneling regime, and a further shortcoming of 
the model is that it cannot account for any variation in 
the distribution with wavelength. It is therefore desirable 
to apply a theory that is expected to be valid also in the 
multiphoton regime. An example of such a theory is the 
molecular strong-field approximation with the inclusion 
of nuclear motion lalSlllIJli and here we apply this theory 
along the lines discussed in 0. This model has previ- 
ously reproduced various experimental observations quite 
successfully 0| . For a linearly polarized laser of fre- 
quency uj and periodicity T = 2tt/lu, the rate of ionization 
differential in the direction q of the ejected electron, to a 
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particular vibrational state of the molecular ion after ab- 
sorption of n photons is dWfi/dq = 2'K\Afi\^qn [atomic 
units, e ~ h ~ = oq = 1, are used throughout]. The 
magnitude of the momentum g„ is determined by energy 
conservation and Afi is the transition amplitude, 

Af, = Sf.^ J\fit)\VF{t)m)dt, (1) 

where Sfi is the FC factor corresponding to the nuclear 
vibrational transition, Vpit) is the length gauge form for 
the molecule-light interaction and \i{t)) and \ f{t)) are the 
initial and final electronic states, respectively, including 
energy phases for the nuclear vibrational motion. The 
time scale of molecular rotations is much slower than the 
typical pulse duration so that the molecular orientation 
is considered as being fixed throughout the pulse. The 
time scale of the vibrational motion on the other hand, 
is shorter than the typical pulse length and therefore we 
treat this degree of freedom quantum mechanically. The 
initial electronic state is the highest occupied molecular 
orbital, and in the final state it is assumed that the laser- 
electron interaction is much stronger than the electron- 
molecule interaction, such that the final state of the elec- 
tron can be accurately described by a Volkov wave. In 
Eq. Q , the electronic matrix element is evaluated at the 
equilibrium distance of the nuclei in accordance with the 
Born-Oppenheimer approximation. This allows us to iso- 
late the overlap between the nuclear wave functions ex- 
pressed by the FC factor. Under experimental conditions, 
the temperature will be sufficiently low so that only the 
vibrational ground state of the initial state is populated. 
We integrate over all directions of the outgoing electron 
and sum over all accessible numbers of photon absorp- 
tions to obtain the total rate of ionization to the vibra- 
tional state considered. Finally, and this turns out to be 
essential, we integrate the total rates over the Gaussian 
spatial and temporal profile of the laser pulse and av- 
erage over molecular orientations for direct comparison 
with experimental data. The spatial integration is car- 
ried out over the interaction region which is restricted 
to a small volume around the beam waist in the experi- 
ment. It is not difficult to show that relative signals are 
independent of focal spot size. 

In Fig. ^ we present the experimental and theoreti- 
cal vibrational distribution of produced by strong- 
field ionization of H2 for two different intensities. For 
comparison also squared FC factors are shown. At both 
intensities, we find good agreement between experiment 
and theory. Both distributions favor the lower vibrational 
states in contrast to the FC distribution. 

Our aim is to provide a widely applicable model and 
therefore we simply use the field-free FC factors. These 
are readily available in the literature for a wide range 
of molecules. In reality, the field-free nuclear potential 
curves are distorted by a very intense laser field, e.g., in 
this effect leads to a bond-softening due to a coupling 
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Figure 1: The vibrational distribution of after ioniza- 
tion in an intense laser field. The black bars indicate the 
experimental observations 4j and the gray bars are the pre- 
dictions according to the present theory. The Franck-Condon 
distribution is indicated with white bars. In both panels the 
laser wavelength is 800 nm, the pulse duration (full width at 
half maximum) is 45 fs. The peak laser intensities are (a) 
3.0 X 10^=* W/cm^ and (b) 4.8 x 10"W/cm^ 



between the lowest Eg and E^, states. As the potential is 
changed, a new set of vibrational eigenstates and corre- 
sponding FC factors needs to be considered. If, however, 
the intensity is below 5 x 10^'^W/cm^ the field-induced 
modification of the potential curve is minimal and the 
application of the field-free FC factors accurate. Also, 
we note that dissociating channels become increasingly 
important as the intensity exceeds 5 x 10"'^^ W/cm^, and 
since one might be interested in making subsequent ex- 
periments on the bound molecular ion we will not consider 
such intensities here. 

Even though the FC principle is explicitly used in our 
theory [see Eq. (Q)], the predicted and observed vibra- 
tional distributions deviate significantly from the FC dis- 
tribution. The key concept explaining this discrepancy is 
the effect of channel closings in connection with full ac- 
count of the pulse profile. By energy consideration, the 
number of absorbed photons n must fulfil the criterion 
nuj = + Ip +Up with Qn > 0, the ionization po- 
tential to the level v and Up = I /{Aur) the quiver energy 
of a free electron in the laser field of intensity /. The pres- 
ence of Up means that the minimum number of photons 
needed to reach the electronic continuum and a particu- 
lar vibrational state v in the ion increases with intensity. 
This phenomena is referred to as channel closing and the 
effect is illustrated in Fig. 12 where we show the ionization 
rates to the lowest vibrational levels of when ionizing 
molecules aligned parallel to the laser field at a wave- 
length of 800 nm. At an intensity of 1 x IQ-^^ W/cm^ all 
the vibrational states shown can be reached by absorp- 
tion of 11 photons. As the intensity increases the thresh- 
olds shift upwards by Up and at the intensities marked 
by arrows, absorption of 11 photons becomes insufficient 
to reach the vibrational levels indicated. For example, in 
the intensity range of 2.3 — 2.7 x 10^'^W/cm^ one needs 
n > 12 photons to reach the v >\ levels whereas n = 11 
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Figure 2: Ionization rates to the lowest vibrational states of 
H J from H2 molecules aligned parallel to the laser field at the 
wavelength of 800 nm. Note logarithmic scale. 



is sufficient for v — Since the rates for higher or- 
der processes (n > 12) are much lower than the rate for 
n = 1 1 , this explains why the j/ = is favored by a factor 
of 3 over v = {1,2}. Contrary, when the same number of 
photons is needed to reach all v levels, the rate to = 
is generally lower than to = {1:2} due to the smaller 
FC factor of the former. 

The reason for finding a very different relative rela- 
tionship between the signals at the peak intensity of 
3 X lO^^W/cm^, Fig.n](a), and the rates at the same in- 
tensity is a result of taking the pulse shape into account. 
Only in the very center of the Gaussian laser beam the 
intensity reaches the peak intensity. In other regions of 
space the molecules are exposed to a lower intensity and 
excitation to i/ = dominates. 

So far we have seen that the vibrational distribution 
depends markedly on the intensity. Our next purpose is 
to investigate how the properties of the laser pulse can 
be varied to maximize the population in a given vibra- 
tional state. A high degree of population transfer to a 
definite state will be a valuable result as it will allow 
for the possibility of making state-specific experiments 
on the molecular ion. To our knowledge, such type of 
control has not previously been explored in strong-field 
physics. We have chosen to model the pulse by a sim- 
ple variation of the peak intensity with the purpose of 
maximizing the v = population. We have performed 
the optimization at the wavelengths corresponding to the 
fundamental- and frequency doubled wavelengths of the 
Ti:Sapphire (400 nm and 800 nm) and Nd:YAG (532 nm 
and 1064 nm) lasers. During optimization the pulse du- 
ration is fixed at 45 fs and 6 ns for the Ti:Sapphire and 
Nd:YAG wavelengths, respectively. 

The result of the intensity optimization is shown in 
Fig. O At the shortest wavelength, 400 nm, a relative 
population as large as 75% is produced in the vibrational 
ground state at an intensity of 2 x 10"'^^ W/cm^. At lower 
intensities we can get nearly 100% in the v = state but 
at such low intensities that the ion yields are very low. In 
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Figure 3: The vibrational distribution of after ioniza- 
tion in an intense laser field according to the present the- 
ory (gray bars). The Franck-Condon distribution is indi- 
cated with white bars. The peak laser intensity is cho- 
sen to maximize the population in the lowest vibrational 
state. The laser wavelengths are (a) 400 nm, (b) 532 nm, (c) 
800 nm, and (d) 1064 nm, and the peak intensities are (a) 
2.0x10^^ W/cm^ (b) 2.7xl0" W/cm^ (c) 2.6xl0" W/cm^ 
and (d) 1.9 x 10^^ W/cm^. The pulse duration is 45 fs in pan- 
els (a) and (c), and 6 ns in panels (b) and (d). 



our model, the reason for obtaining such a confined dis- 
tribution is that the ion yield is completely dominated by 
5-photonabsorption which is only possible to the u = Q 
state and only at intensities below 4.9 x IG^^W/cm^. We 
note that this explanation of the favoring of the v — 
level is different from a previous suggestion offering an ex- 
planation in terms of resonance enhanced ionization [^. 
As the wavelength increases the selection of the = 
becomes less efficient. At all wavelengths the optimum 
intensity is relatively low and it will always correspond to 
an intensity which is slightly below the occurrence of the 
first channel closing. The decreasing v — population 
with increasing wavelength was also observed experimen- 
tally |i,l6(, and the general phenomenon that the largest 
population in the vibrational ground state is obtained at 
relatively low intensities is in good agreement with the 
experiments. 

We note that the control obtained by changing inten- 
sity and/or wavelength is not only the simplest way of 
gaining control, it is probably also the only one. If, 
e.g., one considers the possibility of coherent control by 
preparing either a coherent superposition of initial vibra- 
tional states separated by the energy of one photon or by 
utilizing a bichromatic field and optimizing the strengths 
and relative phase of the two components of the field, one 
faces the problem of a continuum of final states. With 
a proper choice of amplitudes and phase it is possible to 



4 







■ m 


T 


(b) - 


: p 
1 


— 

_i_rn 


- (d) - 

■ yn J1 ^ ,^ , ■ 



1 12 3 4 

Vibrational level Vibrational level 



Figure 4: The vibrational distribution of [(a) and (c)] 
and O2 [(b) and (d)] after ionization in an intense laser field 
according to the present theory (gray bars). The Franck- 
Condon distributions are indicated with white bars. In all 
panels the laser wavelength is 800 nm and the pulse duration 
is 45 fs. A typical peak intensity of 5 x 10^^ W/crn^ is chosen 
in panels (a) and (b). Maximization of the population in the 
lowest vibrational state is obtained with the peak intensities 
of (c) 2.4 X 10" W/cm^ for NJ, and (d) 3.6 x 10^^ W/cm^ for 
0+ 

exclude the transition to one particular final state which 
is characterized by the vibrational state and the direc- 
tion and energy of the outgoing electron. However, one 
cannot gain control over the transition to all other final 
states and hence an efficient selection cannot be made. 

The vibrational wave functions of H2 and lead to 
a quite broad FC distribution and accordingly many vi- 
brational states are populated after ionization. If one 
considers molecules where only a few FC factors are im- 
portant, one may hope for a more efficient optimization. 
To this end, and to illustrate the wide applicability of 
the present formalism, we investigated ionization of N2 
and O2 where the number of final vibrational states are 
limited to < 1 and 4, respectively, simply because the 
other FC factors are vanishingly small. First, we calcu- 
lated the vibrational distribution with a laser wavelength 
of 800 nm, pulse duration of 45 fs and with a typical, but 
arbitrarily chosen, peak intensity of 5 x 10^'^ W/cm^. We 
then made an optimization similar to that for H2 and the 
results are given in Fig. ^ The vibrational distributions 
obtained with a typical peak intensity. Figs. 01 (a) and (b), 
are quite similar to the FC distributions. Figure 0] (a) 
shows typical and Fig.2|(c) optimized vibrational distri- 
bution of after ionization of N2. For N2, only two 
states will be populated in and we see that the FC 
factors strongly favour the = Q state (90%). By utiliz- 
ing a non-optimized strong laser field we obtain a = 
population of 94% [Fig.2](a)], and if we use an optimized 
laser intensity we may reach a. v = population of 97% 
[Fig. 21(c)]. In Figs.0](b) and (d) we present the results 
for O2 and we see that an efficient selection of the 1/ = 
state can also be obtained for this molecule. We find the 



optimal population to be 80% [Fig.2|(d)] compared with 
the FC distribution of 24% and the non-optimized result 
of 30% [Fig. ^ (b)] . When we compare the results of N2 
and O2 with the results of H2 at the same wavelength 
[Fig. El (c)], we see that the former molecules can indeed 
be brought to the i/ = state more efficiently due to the 
limited number of final vibrational states. 

In conclusion, the molecular strong-field approxi- 
mation with the inclusion of nuclear vibrational mo- 
tion explains the recently observed vibrational distribu- 
tions of produced by strong-field ionization of II2 
molecules The theoretical and experimental distri- 
butions are both very different from distributions pre- 
dicted by Franck-Condon factors - a result that can be 
explained by the effects of channel closings in association 
with the spatial and temporal extend of the laser pulse. 
The theory is simple to evaluate and therefore readily 
applicable to diatomic molecules. We proposed a read- 
ily available method of optimizing the population in the 
1^ = state by varying the peak laser intensity. With 
this type of optimization we showed that the vibrational 
distributions vary significantly with the laser wavelength 
in accordance with experimental findings. Finally, we ap- 
plied the theory to N2 and O2 and found that a higher 
degree of control can be obtained than in H2 due to the 
FC factors which are non- vanishing for fewer states in the 
former cases. We note that efficient control of population 
distributions will be of great value in the development of 
state-specific experiments on molecular ions. 
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